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THE DENSITY PROFILES OF MASSIVE, RELAXED GALAXY CLUSTERS: 
II. SEPARATING LUMINOUS AND DARK MATTER IN CLUSTER CORES 
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ABSTRACT 

We present stellar and dark matter (DM) density profiles for a sample of 7 massive, relaxed galaxy 
clusters derived from strong and weak gravitational lensing and resolved stellar kinematic observa- 
tions within the centrally- located brightest cluster galaxies (BCGs). In Paper I of the series, we 
demonstrated that the total density profile derived from these data, which span 3 decades in radius, 
is consistent with numerical DM-only simulations at radii > 5 — 10 kpc, despite the significant con- 
tribution of stellar material in the core. Here we decompose the inner mass profiles of these clusters 
into stellar and dark components. Parametrizing the DM density profile as a power law pdm oc r _/3 
on small scales, we find a mean slope (f3) — 0.50 ± 0.10 (random) fo'13 (systematic). Alternatively, 

cored Navarro-Frenk- White (NFW) profiles with (logr coro /kpc) = 1.14±0.13*q'22 provide an equally 
good description. These density profiles are significantly shallower than canonical NFW models at 
radii < 30 kpc, comparable to the effective radii of the BCGs. The inner DM profile is correlated with 
the distribution of stars in the BCG, demonstrating a close connection between the inner halo and 
the assembly of stars in the central galaxy. The stellar mass-to-light ratio inferred from lensing and 
stellar dynamics is consistent with that inferred using stellar population synthesis models if a Salpeter 
initial mass function is adopted. We compare these results to theories describing the interaction be- 
tween baryons and DM in cluster cores, including prescriptions for adiabatic contraction, and discuss 
possible signatures of alternative DM candidates. 

Subject headings: dark matter — galaxies: elliptical and lenticular, cD — gravitational lensing: strong 
— gravitational lensing: weak — stars: kinematics and dynamics 



1. INTRODUCTION 



The internal structure of dark matter (DM) halos is a 
key prediction of the cold dark matter (CDM) paradigm. 
Numerical simulations follow ing the detailed structure of 
collisionless CDM halos (e.g.,|Navarro et al.|1996|[Ghigna 



et al.| |2000| |Diemand et aj||2nn5t |Graham et al.l|200ti 



Gao et al. |2012bp genericaliy produce a central density 



cusp with pdm ~ t i characteristic of the Navarro- 
Frenk- White (NFW, |Navarro et"aL"1|1996| ) form, proba- 
bly becoming slightl y shallower on very small scales (e.g., 
Navarro et al.|[2010 ). On the hand, simulations are only 
beginning to make predictions for DM halos that include 
baryons, which could profoundly reshape their host ha- 
los. The structure of real DM halos thus contains im- 
portant information about galaxy formation, but there 
is currently no theoretical consensus on the magnitude 
or even sign of these baryonic effects, particularly over 
a wide range in mass. Additionally, the microphysics 
of the unknown DM particle could become important in 
the densest regions, and the inner structure of DM halos 
may there fore provide valuable indirect clues to its na- 
ture ( e.g., |Spergel fe Steinhardt || 2000| lAbaza.jian et al. 



2001 1 |Kaplinghat 2005 ; jPeter et alTpOlDI 

Given the current uncertainty, observations are clearly 
in an good position to guide theoretical efforts. How- 
ever, measurements of DM mass profiles are extremely 
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challenging and are usually limited by confusion with 
baryons, the small dynamic range of the observations, 
and degeneracies that are inherent to individual mass 
probes (e.g., velocity anisotropy). Clusters of galaxies 
are promising locations to make progress. Accurate mass 
measures are available via many independent observa- 
tional probes , especially gravitational lensing and X-ray 
emission (see |Allen et al.|2011| |Kneib fc Natarajan||2011 
for reviews, and references inPaper 1). As we have shown 



(Newman et al. 2009, 2011), combining stellar kinematics 
with strong and weak gravitational lensing yields con- 
straints over 3 decades in radius. This is comparable 
to the best simulations and is thus suitable for detailed 
comparison of the DM profile shape if the baryonic mass 
can be constrained. On small scales in relaxed clusters, 
the latter is dominated by stars in the central brightest 

cl uster galaxy (BCG ). 

Sand et ah] ( |2002| |2004[ ) demonstrated the utility of 



combining resolved stellar kinematics with strong lens 
ing to constrain two-component mass models, i.e., the 
BCG and DM halo separa t ely (s ec Miralda-Es cude l l995" 

~a. 



Natarajan fc Kneib]|1996 ). Sand e"t~aT. l 
clusters and interred a mean {/j) = 0.52 ± 



studied 
.05, where 



-P 



Pdm oc r K , significantly shallower than an NFW cusp 
having /? = 1. They furth er noted pos s ible v ariation in 



(3 from cluster to cluster. Sand et aLl (|2008|) improved 

mi 



on this analysis for two clusters (MS2137 and A383) by 
relaxing the assumption of axial symmetry in the lens- 
ing analysis, instead conducting a full two-dimensional 
study. They found this did not alter their earlier find- 
ings, but noted that the inferred DM slope (3 is sensi- 
tive to the adopted scale radius r s , which could only be 
constrained by additional mass probes at larger radii. 
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Newman et al. 



This was implemented by Newman et al. 
throug h the addition of weak tensing da 



(20091 in A611 



et al. 



In 



Newman 



(2011), we further extended the methodology in 
A383 by constraining the role of pro jection effects (i.e., 
line-of-sight ellipticity ; |Gavazzi|2005[ ) via the comparison 
of X-ray and lensing data. We also presented a radially- 
extended velocity dispersion profile measured in a very 
deep spectroscopic observation. In both A611 and A383, 
we confirmed a shallow inner DM cusp with (5 < 0.3 (68% 
CL) and (3 = 0.59±g jg, respectively^ 

In Paper I of the present series, we presented strong 
and weak lensing and stellar kinematic data for a sam- 
ple of 7 massive (M 200 = 0.4 - 2 x 1O 15 M ), relaxed 
galaxy clusters at z = 0.19 — 0.31. This built upon 
our earlier papers by enlarging the sample of clusters 
with the highest-quality data: weak lensing measured us- 
ing deep multi-color imaging, primarily from the Subaru 
telescope, extended stellar kinematic profiles in the BCG 
obtained primarily at the Keck telescopes, and multiply- 
imaged sources located in Hubble Space Telescope (HST) 
imaging (25 strongly-lensed sources in total, of which 
21 have spectroscopic redshifts). We showed that the 
total inner density profile is remarkably well-described 
by numerical simulations containing only CDM at radii 
> 5 — 10 kpc, despite the significant contribution of stel- 
lar mass on these scales. 

Here we extend Paper I by dissecting the stellar and 
DM contributions, using improved versions of the tech- 
niques developed in our earlier papers. We first show how 
the mass content of the BCGs can be constrained using 
information from the entire sample. We then isolate the 
DM density profiles and quantify their behavior on small 
scales. We show that the DM profiles become shallower 
than NFW models within « 30 kpc, roughly the typi- 
cal effective radius of the BCGs. Furthermore, the inner 
DM density profiles exhibit likely variation from cluster 
to cluster, and this variation is correlated with the prop- 
erties of the BCG. Finally, we interpret our results in 
the context of the recent theoretical literature, focusing 
on the interactions between baryons and DM in galaxy 
clusters and the likelihood that cores in galaxy clusters 
are imprints of DM particle physics. 

Throughout we adopt a ACDM cosmology with f2 m = 
0.3, ft A = 0.7, and H = 70 km s _1 Mpc -1 . Error bars 
and upper limits encompass the 68% confidence interval. 
When pairs of errors are quoted, they refer to the random 
and systematic components, respectively. 

2. DATA AND MODELING 

Whereas the total density profiles were studied in Pa- 
per I, the goal of this paper is to use our two-component 
fits to separate the stellar and dark mass contributions 
in the cluster cores. All aspects of the data and modeling 
were discussed extensively in Paper I (Section 7). Here 
we provide a summary of the features most relevant for 
this paper. Firstly, the stellar mass in the BCG is mod- 
eled based on fits to the surface luminosity in HST imag- 
ing. A uniform stellar mass-to-light ratio T* is assumed 
within each BCG. As discussed in Paper I (Sections 5.1 
and 9.3), this is justified by the mild or null color gradi- 
ents observed over the relevant radial interval. Non-BCG 

4 In the present paper we p rese nt a significantly revised mea- 
surement for A611; see Section 15.31 



cluster galaxies - relevant as perturbations in the strong 
lens model - are included via scaling relations based on 
the fundamental plane (Paper I, Section 7). 

Secondly, the the cluster-scale smooth DM halo is 
parametrized using either a generalized NFW (gNFW) 
model 

Pdm(t) = s (1) 

(r/r s )P(l + r/r s ) 6 P 

or a cored NFW (cNFW) model with 



pDM(r) 



{l + br/r s ){l + r/r s y 



(2) 



Both models have the same large scale behavior (pdm 
r~ 3 ), but the gNFW model contains a central power- 
law cusp with dlogpuM/dlogr — > — /3 as r — > 0, while 
the cNFW model asymptotes to a constant-density core 
within a characteristic radius r corc = r s /b. Both models 
contain the NFW profile in the limits /3 = 1 and r corc — > 
and therefore allow us to explore deviations from canon- 
ical CDM halos in the central regions. Broad, uninfor- 
mative priors are placed on the halo parameters (p s ,r s , 
and j3 or b) and T + (Paper I, Table 7). 

Based on the close alignment between the optical cen- 
ters of the BCGs and both the X-ray centroids (typically 
separated by ~ 3 kpc, comparable to the measurement 
errors) and the lensing-derived centers of mass, we fix the 
center of the halo to that of the BCG (Paper I, Sections 
2 and 7.3). Furthermore, mass estimates derived from 
lensing agree well with independent X-ray observations, 
which constrains the line-of-sight ellipticity of the halo to 
be mi ld in all clusters exce pt A383 (see Paper I, Section 
8 and |Newman et al.||2011| . 

We do not specifically distinguish the hot gas in the 
intracluster medium (ICM), which is thus implicitly in- 
cluded in the halo in our models. Since the distribution of 
the ICM is similar to that of the halo and com prises only 
a ~ 10% mass fraction (e.g., Allen et al.|2004 ), subtract- 
ing the ICM to isolate the DM has very little effect on 
the slope of the density profile (Arflogp/rifogr < 0.05 ; 



Newman et al.|[2~009 Sommer-Larsen fe LimousinpOlO I 
which is the mam focus ot this paper. 

These models are constrained by three data sets. 
Firstly, the mass on scales of ~ 100 kpc to 3 Mpc 
is constrained using gravitational shear (weak lensing) 
measured in deep, multi-color images primarily from 
the Subaru telescope. Secondly, the angular positions 
and redshifts of background galaxies that are strongly 
lensed by the clusters precisely constrain the mass from 
~ 20 — 100 kpc, varying from cluster to cluster. In total 
we located 25 multiply-imaged sources, of which 21 have 
spectroscopic redshifts (7 were first presented in Paper I). 
Finally, the most unique aspect of our analysis is the in- 
clusion of spatially-resolved stellar kinematics within the 
BCGs. These measures are derived from long-slit spec- 
tra primarily obtained at the Keck telescopes. The stellar 
kinematic data typically extend to R rj 10 — 20 kpc and 
display a very homogenous shape that rises with radius, 
indicating a rising total mass-to-light ratio as expected 
at the centers of massive clusters. As demonstrated in 
Paper I (Section 9), the mass models provide good fits 
to the full range of data. 

For the purpose of distinguishing dark and stellar mass, 
the most important physical assumptions are that stel- 
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Fig. 1. — Degeneracy between T t y an d the DM inner density profile when the halo is parametrized with a gNFW (left axis, filled 
contours) or cNFW (right axis, lines) model. The 68% and 95% confidence regions are shown. Arrows at the bottom of each panel show 
estimates T^y S derived from SPS fits to the broadband colors of the BCGs, assuming Chabrier and Salpeter IMFs; the distribution of 

these is shown in the bottom right panel. The permitted range in T*v is set proportionally to T^y and therefore varies slightly from 
cluster to cluster. 

lar mass follows light, as justified above, and that the 
DM halo is adequately described by a gNFW- or cNFW- 
like profile. The precise parametric form is not as criti- 
cal as the presumption that the DM density turns over 
smoothly at small radii - either to a power-law cusp in 
the gNFW case, or to a constant density in the cNFW 
models - without a sharp upturn on small scales. This is 
reasonable: by design these profiles describe pure CDM 
halos in the appropriate limits, and although the effects 
of adding bar yons are uncertain, adia batic contraction 
prescriptions ( Gnedin et al.|2004 2011 ) predict DM pro- 
files that are well- fit by giN b' W models when applied 
to halos and BCGs representative of our sample. Due 
to the density and radial extent of observational con- 
straints (extended stellar kinematic profiles, strongly- 
lensed galaxies usually at multiple redshifts, weak lens- 
ing), we emphasize that we are able to consider quite 
general families of DM halos in each cluster. 



3. SEPARATING LUMINOUS AND DARK MASS: THE ROLE 
OF THE STELLAR MASS-TO-LIGHT RATIO 

In individual clusters there is a degeneracy between the 
stellar mass-to-light ratio T„v = M*/L v and the inner 
DM slope. This is illustrated in Figure [IJ which shows 
results for the mass models summarized in Section [2] and 
derived in Paper I (Section 9). This degeneracy is ex- 
pected, since stellar mass in the BCG can be traded 
against DM. Owing to the multiplicity of constraints de- 
scribed above, particularly kinematic measurements at 
small radii in the stellar-dominated regime, the model 
degeneracy is not complete, and each cluster does carry 
information on both T* and /? or b. 

It is already evident in Figure [T] that most of the clus- 
ters in our sample prefer a DM inner slope that is shal- 
lower than an NFW profile (i.e., f3 < 1), consistent with 



increased if additional information regarding is avail- 
able. Indeed, most clusters are consistent with a wide 
range of when viewed in isolation, due to the uncer- 
tainty arising from the degeneracy described above. Fur- 
thermore, the figure suggests a possible variation from 
cluster to cluster in the DM inner slope, but this con- 
clusion may be contingent upon substantial variations in 
as well. We do not have strong a priori expectations 
about the possible variation from cluster to cluster in the 
DM inner slope, particularly recalling the uncertain role 
of baryons in theoretical predictions. There are, however, 
several strong reasons to believe that the true physical 
variation in within our sample is small. 

Firstly, Figure [l] shows estimates of the stellar mass-to- 
light ratio T^y S derived by fitting stellar population syn- 
thesis (SPS) models to the broadband colors of the BCGs 
(see Paper I, Section 5.2). Currently, SPS models cannot 
predict absolute masses more accurately than a factor of 
2, primarily due to the unknown initial mass function 



(IMF), which we discuss further in Section 5.1 On the 
other hand, relative stellar masses are more robust, espe- 
cially within a homogeneous galaxy population. As the 
bottom right panel of Figure [T] demonstrates, the range 
in T^y S withi n our samp le at a fixed IMF is small. As- 
suming a |Chabrier| ( 12003 ) IMF, the median (T% s 



the full range is on 
9%0 



2.2; 

y ESQ — 2.32, and the rms scatter is 



5 Throughout, Ly arltl ""~*v refer to the observed luminosity, 
including any internal reddening from dust within the BCG. If 
we removed the reddening to obtain the intrinsic Ly and T*y of 
the stellar populations, their scatter would increase. (Reddening 
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Fig. 2. — Probability densities for logasPSi which parametrizes 
the stellar mass-to-light ratio T«y relative to SPS models (Equa- 
tion [3|, are shown for each cluster (thin lines) and jointly for the 
entire sample (thick). The dotted curve shows the effective prior, 
composed of the flat prior on log T,y convolved with a Gaussian 
uncertainty on T^y S of ctsps = 0.07 dex. Arrows indicate the 



effect of adopting mildly anisotropic orbits with a 



: ±0.2. 



Secondly, the rms dispersion in their absolute V-band 
luminosities Ly of the BCGs in our sample is only 
0.1 dex. This small variation is consistent with previ- 
ous studies of BCGs as "standard candles" with uniform 
luminosities and colors (e.g 



, ,., Sandage 1972 Postman fc 

r|1995| [CoIiTns fc Mannjra^Eernaxdi et al.||2007[ ) 

ly, the environments ot the BCGs are the same: 



Lauerl 

f inally, the environments ot the BGGs are the same: by 
construction they are all central galaxies in massive clus- 
ters, and their central velocity dispersions are compa- 
rable. It would be very surprising if this uniformity in 
luminosity and T^y S , which are thought to derive from 
a similar assembly history, were the result of a conspir- 
acy that masks larger variations in stellar mass. Instead, 
based on these physical similarities, it is very likely that 
the BCGs in our sample have si mila r stellar masses and 
T*y As we discuss in Section 5.1 this is further sup- 



ported by recent, independent studies. 

With the well-motivated assumption that the BCGs in 
our sample have a similar T + y, we can use the full sam- 
ple of 7 clusters to jointly constrain its value, thereby 

is indicated only in cool core clusters hosting some current star 
formation.) However, the SPS stellar mass estimates, which are 
significant for our analysis, are much more robust. 



improving the precision and robustness of our measure- 
ments of the DM profile. Before embarking on this, we 
consider how to handle the small variations in T^y that 
we do anticipate, despite the overall similarity. The sam- 
ple spans a redshift range z = 0.19 — 0.31, so some mild 
passive evolution is expected. Additionally, the BCGs 
with the lowest T^y S estimates show optical emission 
lines and far-infrared photometry indicative of ongoing 
star formation (although it involves a small fraction of 
the stellar mass; see Paper I). These BCGs resid e in the 
cool core clusters, consistent with earlier studies (JBildfell | 



et al. 1120081 ILoubser et al. 112009 

1 H I j * — I ' 'J . — " — ' — . 



Therefore, a more precise tec 



Sanderson et al.p 009 ) . 
mique is to define I * for 



each cluster relative to the SPS measurement: 
logasps 



logT,v/T^ v s . 



(3) 



We can then use the full cluster sample to constrain 
(logasps)j which parametrizes a common, systematic 
offset from photometrically-deri ved stellar mass-to-light 
ratios. As describe in Section 5.1 the most probable 
source for large systematic offsets is an IMF that differs 
from that assumed in the SPS models: in this case, that 
of Chabrier. However, our analysis does not depend on 
the physical origin of the offset, only that is it common 
among our BCGs. Since the variation in T^y S is small 
compared to the range of T*y explored in our fits (25% 
versus a factor of 5.3), this approach is not radically dif- 
ferent from assuming a common Y^y. However, it im- 
proves on that assumption by making use of SPS models 
to adjust for small differences in T*y arising from age 
and dust, while making no assumption on the validity of 
their absolute mass scale. 

3.1. Constraining the stellar mass scale 

Figure [2] shows the probability distribution for log «sps 
derived in each cluster. The uncertainty in logasps 
arises from two sources: that in the T*y derived from 
dynamics and lensing, and the uncertainty in T^y S aris- 
ing from random photometric errors. In Paper I we esti- 
mated the latter as ctsps = 0.07 dex. Thus, the probabil- 
ity distributions for log asps are derived by broadening 
those for log T*y by a Gaussian with a dispersion of ctsps- 
We have argued that there are strong a priori rea- 
sons to expect that «sps is uniform across our sample 
of BCGs. Using the probability distributions in Fig- 
ure [2] we can ask whether the lensing and kinematic data 
are indeed consistent with this assumption. One way 
to quantify this is to suppose that the true distribution 
of log asps is Gaussian and infer its intrinsic dispersion 
ciog a- The formalism for inferring the probability dis- 
tributio n P(cTi og g) was dis cussed in Paper I (Section 9, 
and see IBolton et al.||2012|) . The preference for non-zero 



intrinsic scatter can 



hen be assessed by 



AP 



^2ln[P(o loga = CTpcak)/P(cioga = 0)], (4) 



where <T pC ak is the location of the maximum of -P(<Ti ogQ ,). 
For a Gaussian distribution, AP is the number of stan- 
dard deviations from the mean. We find AP = 0.85, i.e., 
a < Ict preference for intrinsic scatter. Thus, the lensing 
and kinematic data are consistent with (although they 
alone cannot prove) our assumption that there is little 
intrinsic variation in asps within our sample. 
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With the physically-motivated assumption that asps 
is the same for each BCG, we can constrain its com- 
mon value simply by multiplying the 7 independent 
probability distributions. The results are shown by 
the thick curves in Figure [2j Very similar values of 
logasps = 0.28 ± 0.05 and 0.26 ± 0.05 are derived using 
the gNFW and cNFW models, respectively, demonstrat- 
ing that these results do not strongly depend on the exact 
halo model. Given the closeness of these results, in the 
following analysis we adopt logasps = 0.27 ± 0.05. De- 
spite marginalizing over fairly general parameterizations 
of the DM profile, we are able to obtain informative re- 
sults due to the high density of observational constraints 
and the sample size. 

Taking the median (T^ s ) = 2.2, we find that 
logasps = 0.27 corresponds to T*y = 4.1. In Section [4~3| 
we describe sources of systematic uncertainty leading to 
a final estimate loga SPS = 0.27 ± 0.05±g;^. In Sec- 
we discuss the physical implications of this re- 



TABLE 1 

Parameters Describing the Inner DM Profile 



5.1 



tion 

suit and compare to the recent literature on the stellar 
mass-to-light ratio and IMF in early-type galaxies. 

4. THE INNER DARK MATTER DENSITY PROFILE 

We now turn to the the inner DM density profiles. In 
our earlier papers, we studied the inner DM density slope 
/3 by marginalizing over the uncertainty in T* separately 
in each cluster. With the benefit of a larger sample with 
improved data, we have now combined constraints from 
7 clusters to arrive at a joint measurement of the stellar 
mass scale asps (Section 3.1 1. Incorporating this infor- 
mation, we can now conduct our analysis in a more phys- 
ically consistent way that recognizes the homogeneity of 
the BCGs, as well as further reducing the remaining de- 
generacies between dark and stellar mass. 

Technically, we implement the joint constraint on 
log agps via importance sampling (e.g., Lewis & Bri- 
dle 2002), reweighting the Markov chain samples to ei- 
iecuvely convert our flat prior on log agps to a Gaus- 
sian with mean (logasps) = 0.27 and dispersion a — 
(a\ + cTgpg) 1 / 2 = 0.09. Here a a = 0.05 dex is the uncer- 
tainty in (logasps)i and ctsps = 0.07 dex is the random 
error in T^y S for each BCG. The latter accounts for the 
fact that asps refers to a systematic offset from SPS- 
based mass estimates, but random errors due to photo- 
metric noise remain in each cluster^ 

4.1. Dark and stellar mass profiles 

Figure [3] shows the resulting spherically- averaged den- 
sity profiles for the DM halo, BCG stars, and their sum. 
The results based on gNFW and cNFW models are again 
quite similar, showing that the choice of parametrization 
does not strongly affect the derived density profiles. We 
do not detect an overall preference for one model over 
the other: the ratio of the total Bayesian evidence is 
consistent with unity. 

The black line segment in each panel spans r = (0.003— 
0.03)r2oo, which is the interval over which the total den- 
sity slope 7tot was defined in Paper I. Its slope r" 



-1.13 



IS 



6 This estimate of o-gpg may be conservative, given that the 
dispersion in T^y S measurements among the BCGs is smaller, and 
X 2 /dof < 1 in the SPS model fits. Thus, in practice we are likely 
allowing for some mild intrinsic variation in asps- 
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Note. — Median parameters are shown, obtained after weight- 
ing samples to incorporate our joint constraint on «SPS: as de- 
scribed in the text. Error bars encompass the 16 — 84th p erce ntiles 
and account for random uncertainties only; see Section |4.3| for a 
discussion of systematic errors. Results are shown for individual 
clusters (top) and for the ensemble mean (bottom), inc luding for 
several alternative assumptions described in Section |4.3| 

the average measured in CDM-only cluster simulations 
from the Phoenix project ( Gao et aL]|2012b ). As quanti- 
fied in Paper I, the stars and DM sum to produce a slope 
very close to CDM-only simulations over this interval. 

Now we can see that both stars and DM contribute 
significantly to the mass in this regime: stars dominate 
the density in the inner radius, while virtually all the 
mass is DM at the outer radius. This demonstrates a 
tight coordination between the inner DM profile and the 
distribution of stars: the NFW-like density slope is not 
a property of the DM halo or the BCG alone, but of 
their sum. As noted in Paper I, at yet smaller radii 
r < 5— 10 kpc where stars are dominant - well within the 
mean effective radius (R e ) = 30 kpc - the total density 
profile generally steepens. 

As expected if the total density is NFW-like, the DM 
profiles become shallower only on scales where the BCG 
contributes significantly, roughly within R e . As we de- 
scribe in Section [5j our results thus do not conflict with 
other studies that claim the DM alone follows an NFW 
profile but are confined to r > R e . The stellar mass den- 
sity in our models reaches that of the DM at a median 
radius of (r) = 7 kpc. In terms of enclosed mass, equality 
occurs at (r) = 12 kpc. Within 5 kpc the median DM 
fraction is (/dm) = 25%, simil ar to massive field ellip- 



ticals (e.g., Auger et al. 2010a), but within their three- 
dimensional halt-light radii r^ the BCGs are far more 
DM-dominated: (/dm) = 80%. 

4.2. Inner DM density slopes and core radii 

Figure [4] shows the probability distributions for j3 
(gNFW) and r core (cNFW) obtained by marginalizing 
over the other parameters, again weighting the samples 
to incorporate our joint constraint on asps- Results for 
the individual clusters are listed in Table[T] Every cluster 
prefers /3 < 1, i.e., an inner slope shallower than an NFW 
model. Thick black lines shows constraints on the mean: 
03) = 0.50 ±0.13 and (logr coro /kpc) = 1.14 ±0.13; the 
method for deriving these is outlined in the Appendix. 



6 Newman et al. 




Radius [kpc] 

Fig. 3. — Radial density profiles for the DM halo, stars in the BCG, and their total in each cluster, adopting the joint constraint 
on log q:imf from Figure [2] The spatial extent of each data set is indicated at the top of each panel. The black li ne segment spans 
r / r 200 = 0.003 — 0.03 and nas the slope p oc r~ 1,13 that is the average of clusters in the Phoenix DM-only simulations I Gao et al.|2012b| 
over the same interval (Paper I, Section 10). The width of the bands indicates the uncertaint y, in cluding both la random uncertainties 
for isotropic models and a systematic component estimated by taking /3 an i so = ±0.2 (Section |4.3| l. Arrows at the bottom of each panel 
indicate three-dimensional half-light radius of the BCG. 




0.4 0.6 0.8 1.0 1.2 
gNFWp (DM inner slope) 




0.5 1.0 
cNFW log r c , 



Fig. 4. — Marginalized posterior probability densities for the gNFW inner slope /3 (left) and the cNFW core radius (right). Thick black 
lines show joint constraints on the mean. Colors follow Figure [2] Each curve is normalized to integrate to unity. Our adopted priors 
(Paper I, Table 7) are flat in logr s and b and therefore not in logr C orc = logr s /b over the whole plotted range. The arrow in the right 
panel indicates the region over which the effective prior on logr coro is flat. 
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y tot : p = -0.75, P = 0.05 
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Effective radius of BCG R e [kpc] 

Fig. 5. — Correlation between the size of the BCG and the inner 
DM profile. Top: Grey points show the total density slope 7 to t 
presented in Paper I; this is measured over r/r2oo = 0.003 — 0.03 
and is not an asymptotic slope. The dashed horizontal line s hows 
the mean slo pe measured in CDM-only cluster simulations | |Gao| 
|et al.|[2012b|l over the same interval. Colored points denote tne 
asymptotic DM density slope /3 measured in the gNFW models. 
Dotted lines show least-squares linear fits. The Spearman rank 
correlation coefficient p and the corresponding two-sided Po -value 
are listed. Bottom: The core radii r core of the cNFW models are 
shown, again indicating a correlation with R e . 

We note that while the typical r core w 14 kpc is small, 
the cNFW profile turns over rather slowly at small radii. 
Thus, while r cor e is the radius where the density falls to 
half of the corresponding NFW profile, significant devi- 
ations extend to r ~ (3 — 4)r corc . 

We can also ask whether there is evidence for intrinsic 
variation in the inner DM profiles. This can be quan- 
tified by assuming that the parent distributions of f3 
and log r core are Gaus sian, and using the method de- 
scribed in Section 3.1 to infer its dispersion. We find 
some evidence for intrinsic scatter with erg = 



and cr\o 



= 0.57 



+0.33 
-0.21- 



n OO+0.15 
u - z ^-0.11 



Its statistical significance can 
be assessed with the AP statistic (EquationEl) : we derive 
AP = 1.5 and 2.6 for /3 and logr coro , respectively. This 
indicates a ~ 2a preference for the presence of intrinsic 
scatter in the inner DM profile shape. 

A possible physical origin of this scatter is illustrated 
in Figure [5j Grey points in the top panel show the total 



density slope 7tot- As described in Paper I, these show 
mild scatter around the mea n slope measured in CDM- 
only simulations (dashed line, Gao et al.|2012b ) over the 
same radial interval (r/r2oo = 0.003 — 0.03). Here we 
see signs of a correlation with the size of the BCG, with 
more extended BCGs corresponding to shallower total 
slopes. The effect on the DM slope (colored points) ap- 
pears stronger: larger BCGs are hosted by clusters with 
shallower DM slopes j3, or equivalently larger core radii 
?"core (bottom panel). Such a correlation is necessary for 
the dark and stellar mass to combine to a similar total 
density profile. The significance can be assessed using 
the Spearman rank correlation test. We find a probabil- 
ities Pn = 0.18 and 0.07 of obtaining an equally strong 
correlation between R e and j3 or r corc , respectively, in 
the null hypothesis of uncorrelated data (see caption to 
Figure [5]) . 

Figure [5] shows that the mass profile in the cluster core 
is closely connected to the build-up of stars in the BCG. 
We return to this point in Section [6] and discuss physical 
scenarios that can explain this. Although the correla- 
tions with R e are most convincing, they are not unique: 
we find correlations between f3 or r core and the stellar 
mass or luminosity with nearly equal statistical signifi- 
cance. There is no sign of a correlation with the virial 
mass M 200 (fi = 0.11 and 0.04 for the gNFW and cNFW 
models; see caption to Figure [5]) F] 

We emphasize that it is preferable to compare directly 
to the physical density profiles (Figure when possi- 
ble, rather than only marginalized distributions for /3. 
These results do not imply, for example, that a CDM 
density profile should be modified simply by maintaining 
the same r s and changing /3 = lto/3 = 0.5. Rather, r s 
also shifts in our fits such that significant changes in pum 
are kept within r < 30 kpc. This degeneracy is simply a 
result of the gNFW parametrization. 

4.3. Systematic uncertainties 

A full discussion of the systematic uncertainty affecting 
our analysis was pres ented in Paper I, Section 9.3 (see 
also Sand et al. 2004 1. In the following, we review the 
most important effects and estimate their impact on asps 
and the inner DM halo parameters j3 and b. 

One of the main sources of systematic uncertainty is 
our use of spherical dynamical models based on isotropic 
velocity dispersion tensors. As discussed in Paper I (Sec- 
tion 9.3), this is a good approximation for luminous, non- 
rotating giant ellipt i cals in their central re gions (e.g., 
IGerhard et"al]|2001| |Cappellari et al " 



2007) 



Nonethe- 
less, individual galaxies can exhibit mild anistotropy with 
|Amiso| = |1 — eg /of I ~ 0-2, and the population as a 
whole also may be slightly radially biased. To estimate 
the impact this has on our analysis, we repeated the dy- 
namical analysis taking a constant anisotropy parameter 
Amiso = ±0.2. Arrows in Figure [2] show that individual 
clusters may shift by AlogT* = —0.16 (f3 aniso = +0.2) 
or AlogT* = +0.10 (/3 an i so = —0.2). Since this bias 
may be correlated among the BCGs, we consider these 

7 Interestingly, the reverse seems to hold for Ttot : there is no sign 
of a correlation with the stellar mass or luminosity, but a possible 
correlation with M200 (p = —0.68, Po = 0.09). The latter may 
simply be because the radial range over which 7tot is measured is 
proportional to T'zoq. 
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as systematic uncertainties in the mean: (log asps) = 
0.27 ± 0.05lg i6- We note that the effects of anistotropy 
are l arger here than for studies of field elliptical lenses 
(e.g., Auger et al.|2010b ), since the latter do not resolve 
kinematics well within K e where the impact of anisotropy 
on the line-of-sight velocity dispersion is largest. 

Uncertainties in the orbital distribution have a milder 
effect on the parameters describing inner DM profile. If 
we adopt the same prior in (logasps), taking /3 a niso = 
±0.2 leads to systematic shifts of A(/3) = ±0.13 and 
A(logr corc ) w -0.18 (Table fT}. If we instead shift the 
prior on (log asps) to match the results obtained with the 
corresponding /3 an i so , we find A(/3) = +0.11,-0.02 and 
A(logr coro ) = —0.21, +0.08. Based on these results, we 
estimate systematic uncertainties of A(/3) = ±0.13 and 
A(logr coro ) = —0.2, +0.1 due to the orbital anisotropy. 

We note that the clusters with the lowest inferred asps 
in Figure [2](MS2137 and A611) are those with the high- 
est halo concentration parameters (Paper I, Section 10). 
These clusters have NFW-like total density profiles down 
to unusually small radii, with very weak steeping on 
small scales. In view of the similarity of «sps among 
the other 5 clusters and the agre ement with independent 
results discussed in Section |5.1[ a likely explanation is 
that some of the stellar mass is effectively counted in the 
halo when Y* is allowed to vary freely from cluster to 
cluster. Nevertheless, omitting MS2137 and A611 would 
shift (logasps) by only +0.02. In this respect our results 
are encouragingly robust. This highlights the utility of 
the ensemble of clusters as a robust constraint on Y*. 

Line-of-sight (l.o.s.) ellipticity in the cluster halo can 
complicate the coupling of lensing and dynamical mass 
measurements, since lensing measures the mass con- 
tained in cylinders, while dynamical and X-ray measure- 
ments nearly measure the spherically-averaged mass dis- 
tribution. The close agreement between lensing- and X- 
ray-based mass measurements shows that this is not a 
major effect in our sample; the only exception is A383, 
in which the l.o.s. shape is explicitl y accounted for in ou r 
analysis (Paper I, Section 8.1, and |Newman et al7||2011 ) . 
Specifically, the mean ratio of spherical mass measures 
(Mx/-Mi ens ) = 1.1 at r ~ 60 kpc, the typical Einstein 
radius in our sample (Paper I, Section 8). This could be 
explained by a mean elongation of the cluster halos along 
the l.o.s. with ellipticity (q — l) ~ 0.1 — 0.2 (although, as 
described in Paper I, (Mx/M\ ens ) and thus (q) are actu- 
ally consistent with unity within the systematic uncer- 
tainties). Based on our study of A383, we estimate that a 
mean l.o.s. ellipticity of this magnitude would cause sys- 
tematic shifts of A(/3) s» 0.06 and A(logr corc ) » —0.1. 

Combining the effects of l.o.s. ellipticity and orbital 
anisotropy in quadrature, we arrive at final measure- 
ments (/3) = 0.50 ± 0.13±g;i| and (logr coro ) = 1.14 ± 
0.13lg22 including random and systematic error esti- 
mates. Naturally, variations in orbital anisotropy or 
l.o.s. ellipticity could cause larger shifts on a cluster-by- 
cluster basis. Such effects could decrease the intrinsic 
scatter in j3 and r core that we infer, but they would have 
to be correlated with the size or mass of the BCG (Fig- 
ure [§. While we have argued that our method of de- 
riving a common value of asps is superior, we note that 
marginalizing over Y*y separately in each cluster as in 
our earlier papers would shift the mean (/?) by < la (see 



, ,, [e = 0.59^37! which does not 



"Separate asps" Table [T]). 

Finally, we recall evidence presented in Paper I that 
A2537 is a possible l.o.s. merger. Such an alignment 
could produce a spuriously shallow DM profile in a lens- 
ing analysis, and A2537 indeed has the shallowest slope 
in our sample. However, Figure [5] provides another expla- 
nation: A2537 has the second-largest BCG in the sam- 
ple. Thus, it does not appear that our results for A2537 
are exceptional. Nevertheless, recognizing its unique na- 
ture in our sample, we note that excluding A2537 yields 

08) = 0.69lg;12 and log r, 
change our main conclusions. 

5. COMPARISON TO PREVIOUS RESULTS 

5.1. Stellar mass-to-light ratio 

These results on the inner DM profile are informed by 
the common stellar mass normalization that we infer, so 
it is important to compare this result t o other measure- 
ments to assess its reliability (see also Cappellari et al. 
2012b for a recent review). As shown iii Section |3.l[ we 
find log asps = 0.27±0.05 for isotropic orbits, witha cor- 
responding T^v = 4.1 ± 0.5 and Y*b = 5.3 ± 0.6 at the 
median T^y S and T^g S . When comparing mass-to- light 
ratios at different redshifts, it is essential to account for 
luminosity evolution. Where necessary, we ev olve sam- 
ples as d\ogT* v /dz = -0.64 ( |Treu et al.|2001| ). We note 
that the ~ 0.05 dex systematic uncertainty in Tfy S (Pa- 
per I, Section 5.2) is relevant only for the interpretation 
of T*v in terms of stellar populations, but it does not af- 
fect the stellar mass and so has no effect on the derived 
mass profiles. 

Discussion of T* is often tied to the IMF. This is be- 
cause the unknown IMF is the dominant source of un- 
certainty in the absolute mass scale for stellar popula- 
tion synthesis models, especially for old galaxies (e.g., 



Bell fc de Jong|2001||Bundy et al.|200 5; Cappellar i et al. 



2006 



|Auger et al.||2009[ |Grillo et al.j | 2009| |Stott et al 
2010 1 . If interpreted as a difference in IMF, our measured 
Qsps ind i cates a normalization consistent with that of the 
Salpeter ( |1955[ ) IMF, which has logM^saip/M^chab — 
0.25 when extended down to 0.1 Mq. 

Several other studies have used lensing and stellar 
dynamics to probe massive field and group ellipticals. 
Auger et al. (2010b) probe the SLACS samples of early- 
type lenses using strong and weak lensing and stella r 
kinematics (see also Gavazzi et al. 2007 Treu et al. 20101). 
Assuming an NFW halo, they inter logasps = 0.28 ± 
0.03 at M* = lO 1:L M p] Assuming an adiabatically- 
contracted halo lowers this value by 0.11 — 0.14, i.e., 
still heavier than a Chabrier IMF. They infer an intrin- 
sic scatter of < 0.09 dex in logaspR w ithin t heir sample 
of cr > 200 km s" 1 lenses dTreu et al.||2010|. iLagattuta 



et al. (20101 study ellipticals at slightly higher redshift 
using strong and weak lensing. Evolving their from 
(z) « 0.6 to our (z) = 0.25 yields Y* v = 4.7 ± 0.7, 
consistent with our results. Both of these works assume 
an NFW halo and a mass-concentration relation that 
follows theoretical expectations (i.e., a one-parameter 
halo). Our models include much more general halos, 
and the BCGs are much more DM-dominated. Thus, 



8 Their oimf is defined relative to a Salpeter IMF and so differs 
from our definition by 0.25 dex. 
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the uncertainty in Y* on an object-by-object basis is 
larger; nonetheless, the ensemble averages agree well. 



and /3dm in the following. 
Most observational studies have focused on the total 



Sonnenfeld et al. (2012) studied a rare early-type lens 
that presents two Einstein rings, which allowed them 
also to relax assumptions on the DM profile. They find 



asps = 0- 3 ± 0.09 in our notat i on (s ee also Spiniello 
et al.||2011[ ) . |Zitrin fc Broadhurstj ( [20091) took advantage 



of the unusually hat surface density profile in the lensing 
cluster MACS J1149.5+2223 (z = 0.544), which offers 
a clean subtraction of the dark halo to isolate the mass 
of the BCG. They estimate T, B ~ 4.5 ± 1 (« 7 ± 2 if 
evolved to our (z) — 0.25). 

Other studies have used intregral field spectroscopy 
to co nstruct detailed dynamical models of local ellipti- 
cals. ICappellari et al.| ( |2012a|b|c| discuss the ATLAS 3D 
sample of early-type galaxies. At the highest velocity 
dispersions prese nt, they infer logasps = 0.25 ( |Cappel-| 
|lari et aL| |2012bj Figure 9, converted to our definition 



of asps)- Interestingly, there appears to be little or no 
intrinsic scatter in asps a t o~ e > 250 km s , nearly at 
the lower limit of our BCGs, although only a handful of 
such objects are present in their sample. Alo ng with the 
tightness of the M/L — a e relation at high a e (Cappellari 



et al.|2012c I , this supports our claim that asp s should be 
nearly const ant within our sample of BCGs. |McConnell| 
et al.'j ( |2011[ ) studied the BCG of A2162 using long-slit 
kinematics and integral field spectroscopy with adaptive 
optics, finding T*r, = 4.6lg'y in their "maximum halo" 
solution. For comparison, our result evolved to z = is 
T* R = 4.1±0.5. 

Finally, the IMF in early-type galaxies has recently 
been studied using detailed spectral synthesis models 
that take advantage of surface gravity-sensitive stellar 
absorption lines. In very high-quality spectra, these con- 
strain the abundance of low-mass dwarfs that contribute 
much to the stellar mass but very little to the integrated 
light. Although the degree of scatter remains unclear, 
these studies suggest that a Salpeter-like IMF - or pos- 



cals (Ivan Dokkum & Conroy 12010 


2012 Conroy & van 


Dokkum 2012 


Smith et ai.||2012). 



In summary, our measurements are consistent 
a variety of other recent works indicating a heavy 
(Salpeter-like) in massive early-type galaxies. En- 
couragingly, studies based on completely independent 
techniques are beginning to converge on the same results. 

5.2. The total inner density slope 

When comparing results on the inner density profiles 
of clusters, it is essential to understand the radial range 
that is being fit and whether the total density profile 
or that of the dark matter is being considered. This 
distinction is most important at radii < 30 kpc where 
the BCG contributes noticeably to the total mass. In 
Paper I we showed that the total density profiles in 
our sample are consistent with CDM-only simulations 
down to r ~ 5 — 10 kpc. The mean total density slope 
(7tot) = 1-16 ± 0.05lg Q7 was precisely measured over 
(0.003 — 0.03)r 2 oo and found to be consistent with colli- 
sionless CDM-only simulations, which have (7tot) = 1.13 
(Paper I, Section 9). Note that 7 to t is measured over 
a specific radial interval and is distinct from the asymp- 
totic inner slopes of gNFW models, which we denote /3 to t 



density profile. Umetsu et al. (2011 ) stacked density pro- 
files for four clusters with high-quality lensing data and 
found that /3 tot = 0.89 ^;^, with the inner 40 kpc//i ex- 
cluded from their fit. |Morandi et al. (2011) measured 
/3 to t = 0.90 ± . 05 in A1689, excluding tlTemner 30 kpc, 



and Coe et al 
tribution is JN 



pOlOl) also found that the total mass dis- 
W-li ke. Using i m aging from the CLASH 



surv ey Q Post man et al] |2012[ ), |Umetsu et~alj ( |2012| ) 
and |Zitrin et al.| ( |2011[ ) derived /3 tot = 0.96+^g (their 
"method 7") andfc~i= 1.08 ± 0.07 in MACS' J1206.2- 
0847 and A383, respectively. These lensing results are 
consistent with our claims that the total density profile 
is NFW-like at r > 5 - 10 kp c. 



Morandi & Limousin ( 2012 1 use lensing and X-ray data 
to derive a total slope /3 tot = 1.02±0.06 in A383 and con- 
trast this with our earlier finding that /3dm = 0.59^35 
in the same cluster ( Newman et al.|201~lj )p| These results 
are not inconsistent. Figure [5] shows that the DM pro- 
file we infer in A383 becomes shallower than an NFW 
model only at r < 30 kpc. These scales are excluded 
by Morandi & Limousin in their fits precisely because of 
the uncertainty in the BCG stellar mass that we have 
addressed using stellar kinematics. At r > 30 kpc the 
total density profile in our models - nearly equal to that 
of the DM - is NFW-like. 



5.3. The dark matter inner density slope 

Among the main scientific goals of studying the inner 
regions of clusters are testing predictions of the collision- 
less CDM paradigm, and understanding the formation 
of the central galaxy and its impact on the DM halo. 
Thus, although precise and robust measurements of the 
total density profile are very valuable, for these goals it is 
clearly important to understand how much of this mass 
is DM and how much is baryonic. Over the past decade, 
we have been developing tools to perform this separation 
fSand et al.|2002||2004[|2008||Newman et al.|2009||2011 b 



21- 



he history of this progress was described in Section [T] 
Sand et al.] ( |2004[ ) measured a mean (/? t ot) = 0.52±0.05 
in a sample oi 6 clusters. We have improved on this ear- 
lier work in many ways: through the use of elliptical lens 
models, the addition of weak lensing data, the incorpora- 
tion of multiple strongly-lensed sources (usually located 
at different redshifts) , the comparison with X-ray results 
to quantify line-of-sight effects, the deeper spectroscopic 
observations of the BCGs that have yielded more precise 
and radially-extended kinematic profiles, and through 
joint constraints on the stellar mass scale asps- This 
work has essentially confirmed our initial findings, with 
the present value (/3dm) = 0.50 ± O.lO+o;^ consistent 



with Sand et al. (2004). (The smaller error bars quoted 
in the latter work are due to the more restrictive model 
assumptions, particularly a fixed scale radius r s .) 

Four of the clusters in the present sample have been 
previously studied in our earlier papers. In gener al our 
result s for MS2137 and A963 are consistent with ISandl 
et al. (2004 2008) within their uncertainties, although 



9 The present measurement of in A383 (Table M jl is slightly 
shallower, but consistent with, Newman et al.| {20lfrdue to our 
new joint constraint on «sps- 
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the present measurements supercede earlier ones due to 
the improvements described above. Our an alysis of A383 
is consistent with |Newman et al.| (2011). The results 



presented here for A6 11, on the other hand, are signif- 



icantly different from Newman et al. 



d2009j): we find 
confidence) . 



0.79l°;^, rather than /3 < 0.3 
This is attributable to two changes in the data: a re- 
vised spectroscopic redshift for a multiply-imaged galaxy, 
and improved stellar kinematic measurements (see Sec- 
tions 4.4 and 6.4 in Paper I). 

As we have shown, it is difficult to separate the BCG 
and DM profiles with lensing alone due to the low den- 
sity (or lack) of constraints near the center. Only in clus- 
ters with exceptional lensing configurations is this feasi- 
ble. An interesting such case is A 1703, whic h presents 
an unusual qua d image close to the BCG. |Limousin| 



et al. ( 2008 1 and Richard et al. ( 2009 1 performed a two- 
component fit - a gJNFW halo and BCG stars following 
light, as in this wo rk - and derive /3dm = 0.92toio4- (See 



Oguri et al 
larger error 



]|2009 

"bar. 



for a consistent result with a much 
This may not be inconsistent with our 



findings, since two clusters in our sample prefer a similar 

see Figure [I} 



slope (A611 and A2390 
be scatter from cluster to cluster 
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and there may 
(120101 



Zitrin et al. 

found that the total density profile 1 in A1703 is well-fit 
by an NFW model. 

X-ray studies of two nearby clusters (A2589 and 
A2029) have also shown that the total de nsity follows 
an NFW profile down to sa 002 - 0.01r vir ( |Lewis et al' 
2003| Zappacosta et al!]|2006 ) . The latter authors noted 
that for any reasonable 1*, this implies a shallower DM 
profile in the central regions where the stellar mass is sig- 
nificant. Their finding agrees well with our work, which 
has quantified the split between stars and DM. [Schmidt 
& Allen (2007) studied a large sample of distant X-ray 
clusters. By assuming a typical BCG stellar mass, they 
estimated (/3 D m) = 0.88±0.29 (95% CL). Often the inner 
~ 40 kpc must be excluded from their analysis, making 
a direct comparison difficult. 

6. DISCUSSION AND CONCLUSIONS 

By combining strong lensing, weak lensing, and stellar 
kinematic observations that extend from ~ 3 kpc to be- 
yond the virial radius, thus spanning the baryon- to DM- 
dominatcd regimes, we constrained flexible, physically- 
motivated models of the dark and stellar mass distribu- 
tions in 7 massive, relaxed galaxy clusters. As discussed 
extensively in Paper I, the density profiles of stars and 
DM sum to produce a slope close to CDM-only simu- 
lations, at least outside the very central « 5 — 10 kpc 
where stars strongly dominate the mass. In this pa- 
per we isolated the dark and stellar density profiles to 
quantify the behavior of the DM on small scales, find- 
ing a mean asymptotic inner power law slope of (/3) = 
0.50 ±0.131 

n.i3i or equivalently a mean DM core radius 



(log 



-0. 

= 1.14±0.13 



+0.14 
-0.22- 



We also presented evidence 



Limousin et al. ] | |2008| imposed a tight prior on the BCG 
stellar mass derived trom Bf S fits, but did not consider uncertainty 
from the IMF. Their SPS estimates are quite high: T^g S 11, 
whereas we find XS| S = 3.0 from fitting the SDSS photometry to 
this BCG, also using a Chabrier IMF. Adjusting the latter to our 
preferred asPS = 0. 27 yields T»r = 5 .7, which agrees with the 
estimate T„b ~ 6 by Zitrin et al.| I 2010^ in this cluster. 
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Fig. 6. — Top: Total density profiles, including baryons and DM, 
for our s ample are overlaid on CDM-only simulations of massive 
clusters ( |Gao et al.||2012b| dashed line, with grey band indicat- 
ing the full range or the simulated clusters; see Paper I, Section 
10). The dot-dashed line shows a system in which an NFW halo 
with concentration C200 = 4.5 is altered using t he modified adi- 
abatic contraction model of |Gnedin et alT] 1 2011[ |. Parameters of 
. An = 1.5, wq = 0.85 were used, with the DOG described by a |Jaffe| 
l |1983[ l profile with scale length rj/r2oo = 0.02 and mass fraction 
M* /M200 = 0.002, which are representative of our sample. The 
radial extent of the data is indicated at the bottom of the panel. 
Bottom: As in the top panel, but showing DM only. (The Phoenix 
simulations thus do not change.) Note that CDM halos match the 
observed total density profiles better than those of DM alone. The 
inclusion of halo contraction (dot-dashed line) only exacerbates the 
difference with the mean observed DM slope (thick black segment). 



for possible variation in the inner DM profile from clus- 
ter to cluster, which correlates with the size and mass 
of the BCG (Figure 5|). This implies a close connection 
between the DM profile in cluster cores and the assembly 
of stars in the BCG. 

The conclusion that the inner DM profile is shallower 
than that of pure CDM halos is fully consistent with 
our previous claims ( S and et al. f2002} |2004l |2008| |New- 



2009 



2011[ ). We have improved upon these 



man et al. 

earlier works by collecting improved data for a larger 
sample of clusters and refin ing our analysis techniques, 
as discussed in Section 5.3 A particular advance en- 
abled by this enlarged, improved sample was a joint 
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constraint on the stellar mass-to-light ratio T* of the 
BCGs in our sample, which we found to be elevated by 

(logasps) = 0.27 ± 0.05t°:i6 dex relative to fits to stel- 
lar population synthesis models that assume a Chabrier 
IMF. Our measurements are instead consistent with a 
Salpeter IMF (or a ny e quivalently "heavy" IMF). As 



reviewed in Section 5.1 this agrees with recent, inde- 
pendent studies of massive, early-type galaxi es based on 
lensing, dynamics, and detailed spectroscopy ( Treu et al. 
20101 [Auger et al.|[20T0bl ICappellari et aL||20l2a|bl |van 
Dokk um fe Conroy||20io1 12012| jUonroy fc van Dokkum 
2012p . These rapid developments in our understanding 
of stellar populations promise significant advances in dis- 
entangling the distributions of dark and baryonic mass 
across a range of systems. 

Figure [6] compares our measurements to high- 
resolution (JDM cluster simulations from the Phoenix 



project (Gao et al. 2012b), clearly demonstrating these 
DM-only simulations are a better match to the total den- 
sity profile than that of the DM alone. In assessing the 
role of baryons on their host halos, much of the focus of 
the theore tical literature has been on the halo co ntrac- 
tion (e.g., |Blumenthal et aL 1986| Gnedin et al.| 



2004 



2011 1 expected to result from a central dissipative build- 
up of baryons. The dot-dashed lines in Figure [6 show 
the effect of ap plying the modified ad iabatic contraction 
(AC) model of |Gnedin et al] pHll to an NFW halo 
and BCG with parameters typical of our sample and of 
the Phoenix simulations (see details in caption). As ex- 
pected, the DM profile steepens (bottom panel), only 
worsening the disagreement with our observations. It 
has been argued that this increase in central DM density 
from AC will boost t he gamma ray flux from DM an- 
nihilation in clusters ( |Ando fc Nagai|[2012[ ) p T l However, 
our results suggest that AC is not the mam process that 
sets the density profile and that the net effect on the 
halo is actually opposite to predictions. (As discussed 
in Paper I, this does not necessarily imply that the same 
theory cannot make valid predictions at the galaxy scale, 
where the star formation efficiency and assembly history 
are very different.) 

The uniformity of the total density slope across a range 
of combinations of DM and stellar density profiles, and 
its close similarity to CDM halos, suggests that the pro- 
cesses that set the inner density profile are mainly gravi- 
tational. A possible formation scenario is that early, dis- 
sipative star formation in the main BCG progenitor cre- 
ates a steep total density slope in the inner ~ 5 — 10 kpc, 
where stars dominate the mass. This size scale is indeed 
similar to the observed sizes of very massive galaxies at 



> 2.5 (e.g., Trujillo et al. 



2006 



van Dokkum et al. 



2008 Newman et al.||20T2 ). Subsequent assembly of the 



extended stellar envelope of the BCG - thought to be 
dominated by low-mass, dry accretion of satellites (e.g., 
Naab et al.l|2009[ |Laporte et~aTl|2012| ) - then mostly re- 
places the DM already in place with satellite material, 
roughly maintaining the density. 

Controlled simulations have indeed shown that dynam- 
ical friction between infalling satellites and the DM halo 
can heat the cusp and reduce the central DM density 



(e.g., |El-Zant et al.|20Qll [2004} |Nipoti et al.|2004[|Jardel 



11 In any case, the highly uncertain contribution from subhalos 
may dominate this si gnal (e.g.,|(jao et al.|2Ul^ap. 



fe Sellwood|2009[|Cole et al.|2011| . This process is dissi- 
pationless, since the orbital energy lost by the satellites 
is transferred to the halo, and thus contrasts with the AC 
picture, in which the baryons' energy is radiated away. 
A connection between the assembled stellar mass and 
the central DM de nsity is naturally expected. Indeed, 



Nipoti et al. (2004) find an anti-correlation between the 
amount of stellar mass assembled in the BCG and the 
inner DM density slope /?, similar to our observations 
(Figure [5j we note the satellites in their simulations in- 
cluded no DM). The strength of the effect depends on 
the d ensity of the satellites and their resistance to strip- 
ping. Laporte et al. (2012) showed that when a stellar 
mass-size relation in line with z > 2 observations is im- 
posed in their simulations (offset by 3 — 5 x in size from 
the local relation), the central DM cusp is flattened to 
(5 ~ 0.3 — 0.7, comparable to our observations. It is im- 
portant to realize that numerical experiments investigat- 
ing this effect have generally lacked a fully realistic and 
consistent treatment of the satellites, so improved simu- 
lations are needed. Nonetheless, the current results are 
promising and likely indi cate that the relevant phys ics is 
primarily dissipationless ( Lackner fe Ostriker||2010 [). We 
note that full hydrodynamicai, cosmological simulations 
that include cooling, star formation, and feedback gener- 
ally do not produce shallow DM cusps or cores (Paper I, 
Section 10), whic h probably reflects k nown overcooling 
effects, although Martizzi et al. (20121 recently showed 
that the inclusion of AGJN feedback greatly improves this 
situation. 

Various DM particle scenarios have also been proposed 
to reduce tension between CDM and observations on 
small scales, including the "missing satellites" problem 
and evidence for ce ntral DM cores or shallow cusps (for 



k|2 

steril e neutrinos at the i ~ key scale (e.g.,[A bazajia n et al 



a recent review, see Primack|2009 ). These include warm 



20011 [Boyarsk y et a.1 |[20q9| |Macci6 et al.||2012| |M 



| |ivienci 

et al.||2012[ ), "f uzzy" CTJMcomposed from an ultr alight 
scalar particle flHu et al.|2000l|Woo fc Chiueh|2009l , dar k 



matter produced from ear fy decays (Kaplinghat 2005 
and dark matter t hat itself decays with a long timescale 
( Peter et al.||2010[ ), among many other possibilties. The 
goal is to preserve the large-scale successes of CDM, while 
allowing for modifications at higher densities where the 
detailed properties of the DM particle might manifest. A 
scenario for which halo density profiles has been worked 
out in detail is a self-interacting dark matter particle 



dSpergel fc Steinhardt||2000l [^shida , et al.||2000 



Dave 



|et al.p00l| . |Rocha et al~1 ( |2Ulij| ) and |Peter et aTf ( |201!j| ) 

showed that a cross-section a ~ 0.1 cm 2 g _i can produce 
« 20 kpc cores in clusters without violating any current 
constraints, e.g., fr om the asphericity o f cluster cores or 



the Bullet Cluster ^Randall et aT|2008| ). Only the dense 
central regions of the halo are affected, where scattering 
can occur within a Hubble time. 

These ~ 20 kpc core sizes are intriguingly similar to our 
observations. On the other hand, they are also very sim- 
ilar to the scale of the baryons, i.e., the size of the BCG. 
It is unclear why the total density profile should then 
match the shape expected of collisionless CDM. In these 
scenarios, the core size arises from the microphysics of the 
DM particle and presumably should not "know" about 
the size of the central galaxy (Figure [5]), for example. 
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Newman ct al. 



Thus, observations of clusters alone cannot provide un- 
ambiguous support for alternative DM theories. Global 
comparisons across a wide range of mass scales (for in- 
stance, a cross-section that also produces correct core 
sizes and densities in dwarf galaxies) remain an essential 
test for attempts to explain low central halo densities in 
terms of the DM particle. 
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APPENDIX 



In Section 4.2 we described how posterior probability distributions P(/3) and -P(log r c 
by weighting the samples in the Markov chains derived in Paper I. The weights 



w 



2na 



exp 



logasps - (loga S Ps) 



3 ) are derived for each cluster 



(1) 



effectively convert a flat prior on logasps (Paper I, Section 7) to a Gaussian with mean (logasPs) = 0.27 and a 
dispersion a = (cr„ +Cg PS ) 1 / 2 . This dispersion accounts for two sources of error: the uncertainty a a = 0.05 dex in the 
global systematic offset (logasps) from SPS estimates Tfy S , and the random photometric uncertainty ctsps = 0.07 dex 
in Tfy S for each cluster. The first uncertainty is correlated across the entire sample, while the second is not. 

Therefore, to obtain constrains on the mean (f3) and (logr colG ), the probability distributions derived for each cluster 
in this manner cannot simply be multiplied, since they are not independent. Instead, we calculate the posterior 
probability of (j3) as 

P«/3))cx y , p((^)|loga S p S )F(logasPs)da S P S . (2) 

Here P{((3)\ logasps) is the posterior distribuion of (/3) at a fixed value of logasps- It is obtained by multiplying 
the probability densities P(f3\ logasps) for the seven clusters in our sample, which are each computed with Gaussian 
weights centered at the fixed value of logasps and a dispersion ctsps (i- e -; a — CT SPS m Equation [T] we now account 
for only the random photometric errors in T^ v s since logasps is fixed). P(logasps), which represents our constraint 
on the common st ellar mass scale, is simply a Gaussian with mean (logasps) = 0.27 and dispersion a a = 0.05 dex, as 
derived in Section [3. 1| for isotropic orbits. 

We estimate the intrinsic scatter in f3 (Section 4.2 1 using the posterior probability densities P(/3|logasps = 0.27) 
for each cluster. That is, we evaluate the cluster-to-cluster scatter in (3 at a fixed value of logasps- All of the above 
comments apply equally to our study of the cNFW models, simply replacing (3 by logr corc . 



